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Abstract 

Two-dimensional (2D) mobile carriers are a wellspring of quantum phenomena. Among vari¬ 
ous 2D-carrier systems, such as field effect transistors and heterostructures, polar materials hold 
a unique potential; the spontaneous electric polarization in the bulk could generate positive and 
negative 2D carriers at the surface. Although several experiments have shown ambipolar carriers 
at the surface of a polar semiconductor BiTel, their origin is yet to be specified. Here we pro¬ 
vide compelling experimental evidences that the ambipolar 2D carriers at the surface of BiTel 
are induced by the spontaneous electric polarization. By imaging electron standing waves with 
spectroscopic imaging scanning tunneling microscopy, we find that positive or negative carriers 
with Rashba-type spin splitting emerge at the surface correspondingly to the polar directions in 
the bulk. The electron densities at the surface are constant independently of those in the bulk, 
corroborating that the 2D carriers are induced by the spontaneous electric polarization. We also 
successfully image that lateral p-n junctions are formed along the boundaries of submicron-scale 
domains with opposite polar directions. Our study presents a novel means to endow non-volatile, 
spin-polarized, and ambipolar 2D carriers as well as, without elaborate fabrication, lateral p-n 
junctions of those carriers at atomically-sharp interfaces. 

PACS numbers: 73.20.At, 68.37.Ef, 71.70.Ej 
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I. INTRODUCTION 


The spontaneous electric polarization of polar materials causes surface charges, or band 
bending near the surface. If the band bending is large enough for the conduction and valence 
bands to cross the Fermi level, two-dimensional (2D) positive carriers can emerge at one side 
of a material and negative ones at the other side. Such mobile carriers, if realized, expand 
an arena of 2D-carrier systems that exhibit a wide variety of quantum phenomena^A 
pioneering work on BaTiOs claimed that the electric conductivity was increased by the 
polarization^. However, carrier polarities were not shown and thus their relationship with 
the polarization is unclear. To clarify the relationship, it is necessary to carefully examine 
surface modifications, as emphasized in Ref. 3. For example, oxygen vacancies created on 
the fractured surface produce a 2D electron gas at the surface of SrTiOs (Ref. 4 and 5). In 
contrast to the simple mechanism, it has never been established whether the spontaneous 
electric polarization actually induces 2D carriers. 

A polar semiconductor BiTel is an emergent candidate possessing the polarization- 
induced 2D carriers. BiTel has a layered crystal structure with triple layers composed of 
Te, Bi, and I layers stacking along the c-axis®"®, as shown schematically in Fig. 1(a). In this 
polar crystal structure, the spin degeneracy in the band structures is lifted by spin-orbit in¬ 
teraction®. Angle-resolved photoemission spectroscopy (ARPES) reveals that BiTel indeed 
exhibits momentum-dependent (Rashba-type) spin splitting larger than ever reported^®“^^. 
Because of this feature, BiTel is proposed as a component for spintronics applications^®. 
In addition, both n-type and p-type two-dimensional (2D) band dispersions are observed 
at the surface of BiTeR^’^^. Although preceding studies have shown the presence of band 
bending near the surface accompanied with the surface carriers^®”^^’^^, the origin of the band 
bending (and thus the origin of the surface carriers) remains elusive. Ishizaka el al. indicate 
similarity to near-surface electron-accumulation layers of semiconductors^®. Eremeev et al. 
suggest that the breaking translational symmetry at the surface with the strong ionicity 
modifies the electrostatic potential near the surface^^. In contrast to these surface origins, 
Butler et al. speculate that the spontaneous electric polarization in the bulk causes the 
spectral shift to discuss the surface structures^®. Actually, the surface structure necessary 
to identify the origin is also still obscure. Multiple termination layers hosting the ambipolar 
carriers are attributed to stacking faults (an excess or deficiency of an atomic layer) and 
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steps^^, or opposite stacking sequences^^’^®’^®. Surface vacancies or absorbates have not been 
addressed in the preceding microscopy 

In this study, we substantiate that the ambipolar 2D carriers at the surface of BiTel 
are induced by the spontaneous electric polarization. A core challenge, defined by the 
previous studies, is to elucidate local electronic states including carrier polarities, together 
with surface structures from submicron- to atomic- scales. For this purpose, we performed 
spectroscopic-imaging scanning tunneling microscopy (SI-STM) which yields images of the 
local density of states (LDOS) by measuring spatial variation of the differential tunneling 
conductance, d//dD. To determine local carrier polarity, we exploit electron standing waves 
appearing as spatial modulations in d//dD images^°’^h An electron standing wave is caused 
by quantum interference between electron waves incident to and elastically scattered from 
an atomic defect or step. The wavevector of electron standing wave is given as difference 
between those of original electron waves, q = kg — kj. Therefore, the dispersion relationship 
q{V), specifically the sign of d|q|/dl/, directly reflects carrier polarity: negative for holes 
and positive for electrons. 


II. METHODS 

Single crystals of pristine and substituted BiTel used in this study were grown by a 
modified Bridgman method®. All samples are doped with electrons due to non-stoichiometry. 
Bulk electron densities are determined by the Hall coefficient measured at room temperature. 
The 10^® cm“® samples are metal with temperature-independent electron densities®’^^. The 
Fermi levels of 10^® cm“® samples lie ~0.1 eV above the bulk conduction band minimum^®“^^. 
Meanwhile, since the 10^^ cm“® samples are not fully metallic and have the Fermi levels 
near the bulk conduction band minimum, their electron densities at low temperatures where 
SI-STM measurements were performed may be smaller than those at room temperature. 
Therefore, the 0.1 eV difference of the Fermi levels and the two orders of magnitude difference 
of the bulk electron densities are the minimum estimates of those at low temperatures. 

For SI-STM measurements, BiTel crystals were cleaved in an ultra-high vacuum chamber 
at ~77 K. The crystals were then immediately transferred with a transfer rod cooled together 
with the crystals, through an insert cooled by liquid helium, to a home-built STM head 
placed at the bottom of the insert and cooled down to 4.6 K beforehand. All SI-STM 
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measurements were carried out at 4.6 K with tungsten tips sharpened electrochemically and 
prepared with a field ion microscope. Bias voltages were applied to the sample and the 
tip was virtually grounded. Topographic images were taken in the constant-current mode. 
d//dl^ spectra and images were measured by a standard lock-in technique with a setup 
current of 0.2 nA and the feedback loop opened. The modulation voltage was 5 mVrms 
unless otherwise noted. Fourier transforms of dl/dV images are symmetrized based on 
symmetry of BiTel. 


III. RESULTS 

To grasp submicron-scale characteristics of BiTel, we show a large topographic image in 
Fig. 1(b). A prominent structure observed in all samples studied is the two types of domains 
identified topographically and electronically as described below. The domains are typically 
several hundred nanometers in size and separated by depression in topographic images. The 
domain boundaries are not straight along the crystalline axes but meandering smoothly. 
One type of the domains is apparently higher than the other by about 0.1 nm, which weakly 
depends on bias voltages (Appendix A) . The bias dependence means that the topmost layers 
of the domains are crystallographically inequivalent. 

A topographic distinction between the two types of domains is defect patterns appearing 
in high-resolution topographic images shown in Fig. 1(c) and 1(d). Besides the triangular 
lattices with the bulk a-axis constant, atomic defects are imaged as three dots in Fig. 1(c), 
but dark triangles in Fig. 1(d). The fact that only two kinds of domains are observed 
indicates that the domains are composed of opposite stacking sequences, I-Bi-Te (I-top) and 
Te-Bi-I (Te-top). (See Appendix B for detailed description.) More details about the domain 
structure are brought by further investigating topographic images. First, stacking sequences 
of each domain are I-top for Fig. 1(c) and Te-top for Fig. 1(d). Second, the structural 
relationship between the domains is 2 (reflection about the (0001) plane). Third, central 
sites of defect patterns are identified as shown in Fig. 1(c) and 1(d). (See Appendix C for 
full descriptions.) These thorough identifications of the domain structure play an important 
role to understand electronic states of each domain as described below. 

An electronic distinction between the two types of domains is represented by dl/dV 
spectra shown in Fig. 1(e). The spectrum of the Te-top domain shows finite conductance 
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Figure 1. (Color online) Two types of domains of BiTel. (a) Crystal structure of BiTel. (b) 
A 0.49 /im X 0.28 /im topographic image. The scanning parameters are 10 pA at -0.2 V. The 
gray arrows with indices show the lattice orientation determined by the Bragg peaks in Fourier 
transforms. Filamentary structures indicated by the other arrows are only occasionally observed 
and irrelevant to our arguments, (c and d) 30 x 30 nm^ topographic images of each domain. The 
scanning parameters are 0.2 nA at -0.2 V for (c), and 0.2 nA at -20 mV for (d). The arrows indicate 
examples of defects. The annotation Dx denote defects at X-site (X = Bi, Te, and I). Dark patches 
in (c) are clusters of some defects, (e) dl/dV spectra averaged in each domain. The setup bias 
voltage is -0.2 V. (f and g) dl/dV images taken in the same areas as (c) and (d), respectively, (f) 
was taken at -195 mV and (g) at -hl95 mV. The setup bias voltages are -0.39 V for (f) and -0.2 V 
for (g). 

in the voltage range studied whereas that of the I-top domain shows vanishingly small 
conductance at positive bias voltages. The latter implies that, at the surface of the I-top 
domain, the top of valence band is located slightly above the Fermi level and thus the charge 
carriers are holes in contrast to electrons in the bulk. dl/dV images also clearly manifest 
the distinction as shown in Fig. 1(f) and 1(g). Electron standing waves are observed mainly 
near the defects on the I-top domain while all over on the Te-top domain^®. By comparing 
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these dl/dV images with the topographic images taken in the same location (Fig. 1(c) and 
1(d)), we find that strong scattering centers are defects at Te-site in the I-top domain and 
Bi-sites in the Te-top domain. 

To unambiguously conclude local carrier polarity of each domain, we focus on the disper¬ 
sion relationships of the electron standing waves by analyzing Fourier transforms of d //dV 
images as a function of bias voltages. In the I-top domain, a branch fans out from the F 
point with decreasing bias voltages as shown in Fig. 2(a) and 2(b)^^. The high-symmetry 
linecuts shown in Fig. 2(c)-2(f) exhibit a branch with negative slope (d|q|/dl^ < 0) and 
crossing 0 mV as well as one more branch in deeper energies. Therefore, as implied by the 
d//dV spectrum shown in Fig. 1(e), a 2D hole gas is formed at the surface of the I-top 
domain. Meanwhile in the Te-top domain, the dispersion relationship is completely different 
from that in the I-top domain. As shown in Fig. 2(g) and 2(h), two dispersive branches are 
observed. The outer branch appears as a peak in the F — M direction and the inner one 
as a hexagonal ring surrounding the F point. Both branches approach to the F point with 
decreasing bias voltages^^. The high-symmetry linecuts shown in Fig. 2(i)-2(l) exhibit that 
these branches have positive slope (d|qf|/dV > 0) and cross 0 mV. Therefore, the surface 
carriers of the Te-top domain are electrons. An energy-independent feature appearing on the 
outside of the outer branch is an extrinsic feature inherent to STSTM because its location 
changes depending on setup bias voltages^®. We note that the identification of the top layer 
and corresponding dispersion relationships are consistent with the ARPES results^^. 

The clear electron standing waves in the Te-top domain bear closer analyses to unveil 
the nature of the electronic state. The presence of two branches indicates that two scat¬ 
tering channels connecting non-orthogonal electronic states are involved with the electron 
standing waves^®”^^. By comparing the observed dispersions with the ARPES results^°“^^, 
we assign the scattering channels in the Te-top domain as illustrated in Eig. 2(m). The in¬ 
ner branch arises from inter-band scattering between the spin-split bands as had ever been 
observed for Rashba-split states^^’^^. The outer branch originates from intra-band scattering 
in the hexagonally-warped outer band. This branch emerges, as known for the surface Dirac 
Fermions of topological insulators^^^^^, because of a non-zero out-of-plane component of the 
spin polarization that is characteristic of (Pav-symmetry of the crystal lattice^^’^®. Given 
this assignment of the two branches, we can estimate the momentum offset of the spin-split 
bands at (^out ~ 9in)/2 0.5 nm“^. This value agrees well with the ARPES results^°“^^. 
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Figure 2. (Color online) Fourier analysis of the electron standing waves, (a, b, g, h) Fourier 
transforms of dl/dV images. The setup bias voltages are -0.2 V for (a) and (b), and 0.2 V for (g) 
and (h). The one-sided gray arrows in (h) and (i) indicate the extrinsic feature described in the 
text, (c, d, i, j) High-symmetry linecuts of the Fourier transforms, taken along F — M (c and i), 
and F — K (d and j). (e, f, k, 1) Peak positions extracted from (c), (d), (i), and (j), respectively, as 
quantitative eye guides to show the dispersive branches. The markers and the error bars are peak 
positions and widths extracted by Lorentzian fits, respectively. The error bars also indicate fitting 
directions. Peak positions extracted from another dataset taken with a setup voltage of -0.2 V are 
also shown in (k) and (1), demonstrating that the dispersive branches are independent of setup bias 
voltages, (m) Schematic figures of the surface band structure of the Te-top domain observed by 
ARPES: from top to bottom, the Fermi surface, a 3D illustration of the band structure, and the 
dispersion along F — M. The gray arrows and markers denote directions in-plane and out-of-plane 
spin components. The latter component arises concomitantly with the hexagonal warping. and 
gfout are scattering vectors of the observed electron standing waves, ko is the momentum offset of 
the spin-split bands. 
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Figure 3. (Color online) n-independence of the 2D electron gas induced by the spontaneous 
electric polarization, (n is bulk electron density.) (a-c) Fourier transforms of dl/dV images taken 
at 0 mV with a lock-in modulation of 1 mVrms and a setup bias voltage of -20 mV. The solid white 
lines enclose the doubled g-space Brillouin zone where scattering inside the first Brillouin zone of 
fc-space is found. The broken orange lines are trajectories along which the linecuts shown in (d) 
were taken, (d) F — M linecuts of the Fourier transforms shown in (a)-(c). Each curve is shifted 
vertically for clarity, gin, ^out? and qo indicate positions of the inner branch, the outer branch, and 
the Bragg peak, respectively. The peaks annotated as qo-Qout and qo-qin are replicas of the outer 
and inner branches, respectively, (e) A band diagram along the c-axis of BiTel and a corresponding 
schematic of the crystal structure. Ec, Ey, and Ep are the bottom of conduction band, the top 
of valence band, and the Fermi level, respectively. The arrows annotated with “P” indicate the 
direction of the spontaneous electric polarization. 

The spin splitting is observable in the electron standing waves of spin-split bands (without 
scattering to additional bands^^’^^) when the band structure deviates from a simple Rashba 
model of a theoretical prediction^^ and is hexagonally warped. 

The most salient feature of the electron standing waves in the Te-top domain is found in 
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relationship to bulk electron density, n. As shown in Fig. 3(a)-3(d), the Fourier transforms 
of dl/dV images reveal that the wavevectors of two branches are independent of n. This 
feature was commonly observed in the samples studied and thus is also independent of details 
of the domain structures. The constant wavevector means that, as the Fermi level shifts as 
n changes (~0.1 eV for the n range studied^^), the surface band also shifts such that the 
Fermi wavelength of the surface 2D state stays constant. As indicated by the constant Fermi 
wavelength, the electron densities at the surface stay constant, even though those in the bulk 
change by two orders of magnitude. The surface electron density is roughly estimated at 
(9in/2)V(27r) ~ 3 X 10^^ cm The peculiar n-independence is essential for specifying the 
mechanism of the 2D carriers as discussed below. 


The electron standing waves are also observed around the domain boundary, providing 
details of the p-n junction at the boundary. Fig. 4(a) and 4(b) show a topographic and a 
dl/dV image taken in the same area around a domain boundary^®. To highlight atomic-scale 
structures, we apply a Fourier filter suppressing long-wavelength structures. A structural 
model of the boundary is suggested from a filtered topographic image shown in Fig. 4(c). 
(See Appendix D for details.) The electron standing waves in both domains are better 
visualized in a filtered dl/dV image shown in Fig. 4(d)^®, indicating that this p-n junction 
is a tunnel diode between highly-doped semiconductors. The depletion layer of the p-n 
junction manifests itself as the intervening zone without electron standing waves in Fig. 4(d). 


Transition of the electronic states across the p-n junction is summarized in Fig. 4(e). 
dl/dV spectral variation occurs in the intervening zone, defining the depletion layer. The 
width of depletion layer is about 6 nm, in agreement with a simple estimate. (See Appendix E 
for details.) The depression observed in the topographic image corresponds to the depletion 
layer. The work function has three characteristics; it has larger value in the Ftop domain, 
shows larger modulations in the I-top domain, and changes solely in the Te-top domain side 
of the depletion layer. The larger value is related to the origin of the ambipolar 2D carriers 
as described later as well as apparent heights of the domains (Appendix F). Further analysis 
is required for the larger modulations and the asymmetric change. 
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Figure 4. (Color online) Electronic structures of a lateral p-n junction of the 2D carriers, (a) A 
20 X 20 nm^ topographic image around a domain boundary. The scanning parameters are 0.2 nA 
at -0.1 V. (b) A dl/dV image taken in the same area as (a) at 0 mV with a setup bias voltage of 
-0.3 V. (c and d) Fourier-filtered images of (a) and (b), respectively. Low-wavenumber components 
are suppressed to enhance contrast of atomic corrugations. The broken lines are a guide to eyes for 
showing lateral shift between the topmost sublattices of the two domains. The insets are zoom-in 
images in the areas of the boxes, (e) Intensity plots of dl/dV (top) and I-z (middle) measurements 
taken along the trajectory shown as solid lines in (a), (b), and (d). A setup bias voltage of -0.3 V 
and a lock-in modulation of 4 mVrms are used for the dl/dV measurement. The bottom panel 
shows the corresponding profile of (a) and the work function calculated from the I-z data. The 
positions 1 and 3 are edges of the depletion layer, and the position 2 is the middle point between 
1 and 3. These positions are indicated as the circles in (a)-(d) and the broken lines in (e). 

IV. DISCUSSION 

How can the ambipolar 2D carriers be induced at the surface of BiTel? The key ingredient 
to specify the origin is the peculiar constant density of the 2D electrons. Similar constant 
density is reported for SrTiOs where the 2D electrons are attributed to surface oxygen 
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vacancies^. In this case, the density of oxygen vacancies at the surface is much larger 
than and virtually independent of that in the bulk. In contrast, however, extrinsic surface 
modifications are not observed in our samples (Appendix G) and therefore excluded as the 
origin of the 2D carriers. Instead, the constant density in BiTel is naturally explained 
by the spontaneous electric polarization in the bulk. The unit BiTel layer consists of a 
positively-charged (BiTe)+ bilayer and a negatively-charged 1“ layer^. The spontaneous 
electric polarization therefore directs from a BiTe bilayer to an I layer within the unit layer. 
The conduction (valence) band then is bent negatively (positively) and split off to form a 
2D electron (hole) gas at the surface of the Te-top (Ttop) domain, as illustrated in Fig. 3(e). 
Since the density of accumulated 2D electrons is determined by strength of the polarization, 
it is independent of bulk electron density. 

The other observations are also comprehensively understood by the spontaneous electric 
polarization. The work function larger in the I-top domain is consistent with stronger con¬ 
finement of electrons into the bulk caused by the spontaneous electric polarization pointing 
from a BiTe bilayer to an I layer. The strong scattering at Bi-site (Te-site) in the Te-top 
(I-top) domain reflects the orbital character of the surface band split from the conduc¬ 
tion (valence) band where Bi-6p (Te-5p) orbitals predominantly contribute^®. The domain 
structure with the 2 relationship and the meandering boundaries is analogous with those 
of displacive-type ferroelectrics. We thus conclude that the ambipolar 2D carriers at BiTel 
surfaces are induced by the spontaneous electric polarization in the bulk. 

In this study, we have precisely determined local carrier polarity by observing electron 
standing waves with STSTM. This results demonstrates that, to probe carrier polarity, 
SI-STM is available complementarily to other techniques such as ARPES, thermoelectric 
probes, single electron transistors, and photoelectron emission microscopy. Moreover, we 
unveiled that the underlying mechanism of electron standing waves is common to that of 
Dirac Fermions in topological insulators. This enables us to examine Rashba-type spin 
splitting with STSTM and provides a unified framework to address spin-dependent scatter¬ 
ing phenomena that are a key aspect for spintronics applications. Most importantly, our 
study establishes that the spontaneous electric polarization induces ambipolar 2D carriers. 
Carrier densities of these 2D carriers may be controlled by strain and temperature via the 
piezoelectric and pyroelectric effects. Since the p-n junction at the domain boundary is a 
consequence of the ambipolar 2D carriers and the domain structure sharing the common 
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Figure 5. (Color online) Bias dependence of apparent height. The orange curves are line profiles 
taken across a boundary of the two types of domains. The background image shows a 20 nm x 
8 nm topographic image taken at -0.1 V around the boundary. The white line denotes a trajectory 
where the profiles were taken. 

root, the surface of a polar material is a new platform to investigate a lateral p-n junction of 
2D carriers. Such a surface can be a potential substructure to study unconventional devices 
such as a topological p-n junction^^. The spontaneous electric polarization is, as a means 
to induce 2D carriers, complementary to field effect transistors and heterostructures in the 
sense that the induced carriers are non-volatile and polarity-switchable, and thus is available 
to explore new phenomena and functionalities. 
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Appendix A: Bias dependence of relative heights 

Apparent heights of the two types of domains depend on the bias voltages as shown in 
Fig. 5, indicating that the electronic states at the surfaces of the domains are different and 
the electronic difference contributes to the difference of apparent heights. This is contrasting 
to a step where its height is independent of the bias voltages. See also Appendix F about 
contributions to the apparent height. 


13 







0 20 40 60 


Lateral distance (nm) 

Figure 6. (Color online) Atomic steps with heights of multiples of the c-axis constant. The solid 
curve is a line profile taken across atomic steps. The background image shows a 74 nm x 31 nm 
topographic image taken at -0.2 V. The dotted line denotes a trajectory where the profile was 
taken. 

Appendix B: Domains with opposite stacking sequences 

The experimental fact that only two kinds of domains are observed means that number of 
surface structures is two. If there are multiple cleavage planes in the ideal crystal structure, 
more than two kinds of surface structures appear due to combinations of the topmost and 
subsurface layers. Stacking faults also give more than two kinds of surface structures for 
the same reason. Multiple cleavage planes and stacking faults are excluded also by the step 
heights that are multiples of the c-axis lattice constant (Fig. 6). Then, the topmost layer 
of each domain must be one of the two layers adjacent to the natural cleavage plane, and 
the second-topmost layer must be one remaining layer. Namely, Te and I layers are the 
top and a Bi layer is the second. That is, the observed domain structures are composed of 
opposite stacking sequences, Te-Bi-I (Te-top) and I-Bi-Te (I-top). This domain structure 
naturally accounts for that spin polarization is observed to be unchanged even when a crystal 
is flipped^^, as suggested by the authors. 


Appendix C: Details about the domain structure 
1. Identification of stacking sequence 

To identify stacking sequences of the domains, we measured BiTei_a;Sea,I where Se sub¬ 
stituted for Te works as a marker of Te site. As shown in Fig. 7(a) and 7(b), Se atoms 
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Figure 7. (Color online) Topographic images of substituted samples to identify the domain 
structure, (a)-(d) 14x14 nm^ topographic images of BiTei-a^Se^I ((a) and (b)) and Bii-a^Aga^Tel 
((c) and (d)). Setup bias voltages are -20 mV (a), -0.4 V ((b) and (d)), and -0.2 V (c). For each 
sample, the images were taken on a single surface. Arrows in (a) indicate several exemplary Se 
atoms substituted for Te atoms. The annotations denote defects at X-site (X = Bi, Te, and 
I), (e)-(k) Close-up images around impurities and defects in 3x3 nm^ squares. Areas of (e)-(h) are 
shown as boxes in (c) and (d). (k) is an image around indicated by the arrow in Fig. 1(c). This 
image is rotated to fit the other images. (1) Point group symmetries to realize opposite stacking 
sequences. Schematic figures of the topmost Te and I atoms are superimposed on topographic 
images of Agsi clipped from (e) and (i). Location of Agsi relative to the topmost sublattice proves 
that the structural relationship between the domains is 2 (reflection). 

are observed at the atomic site of Fig. 7(a), indicating that this type of domain with dark 
triangles has the Te-top stacking and the other with three dots has the I-top stacking. 


2. Point-group operation between the domains 

Here we consider the structural relationship between the two types of domains in terms of 
a point-group operation. As is evident from two topographic images taken on a single surface 
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(Fig. 7(a) and 7(b), and Fig. 7(c) and 7(d)), the orientation of the topmost sublattices is 
identical for the two types of domains. Point-group operations to realize opposite stacking 
sequences with keeping orientation of the topmost sublattices are 1 (inversion), 2 (reflection), 
3, and 6. Actually, only 1 and 2 are independent because of the three-fold crystal symmetry 
of BiTeL 

A distinction between 1 and 2 is given by locations of Bi atoms. With respect to the 
topmost sublattice, there are two possible sites where Bi atoms can occupy. If 1 is the case, 
Bi atoms occupy the same site in both domains. If 2 is the case, they occupy different 
sites. To make this distinction, we measured topographic images of Bii_ 2 ;Ag 2 ;TeI as shown 
in Fig. 7(c) and 7(c). Since Agei atoms are observed at different sites (Fig. 7(e)and 7(i)), the 
point-group operation is 2 (reflection). The above discussion is summarized in Fig. 7(1). This 
successful identification of the domain structure demonstrates that a combination of crystal 
growth and scanning tunneling microscopy works complementarily to diffraction techniques 
and transmission electron microscopy to determine local structures. 


3. Identification of defect patterns 

Fig. 7(e)-7(k) show topographic patterns of defects and impurities. Identifying a defect 
pattern centered at the topmost Te atoms is straightforward (Fig. 7(g)). Although defects 
bright in color are often found at the topmost I site, we have not yet identified a feature 
common to all images of the I-top domain. Bi-site defects (Fig. 7(f) and 7(j)), found at 
an interatomic site, are identified based on locations of AgBi (Fig. 7(e) and 7(i)). Patterns 
centered at the other interatomic site are identified as remaining ones, I-site defects for the 
Te-top domain (Fig. 7), and Te-site defects for the I-top domain (Fig. 7(k)). 


Appendix D: A structural model of the domain boundary 

Supplementary Figure 8 shows a structural model of the domain boundary. This model 
is constructed to satisfy the 2 relationship between the domains and the lateral shift of 
the topmost sublattices shown in Fig. 4(c). As for a possible vertical shift of the domains, 
we can not separate a morphological height difference from the apparent height difference 
shown in the topographic image (Fig. 4(a)) because the LDOS and the work function also 
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Figure 8. (Color online) A structural model of the domain boundary, seen from the top and the 
side. The arrows annotated with “P” indicate directions of the spontaneous electric polarization. 

contribute to the apparent height. (See also Appendix F for a mathematical description.) 
The Bi sublattice, therefore, is drawn as flat for simplicity. 


Appendix E: Width of the depletion layer 


Here we calculate the width of the depletion layer of a planer p-n junction, d, with 
reasonable parameters to compare with the width of the observed depletion layer, d is given 
by 


d = 


l2eoeV£) + N] 


D 


(El) 


V e IVaJVd 

where eo and e are the vacuum permittivity and the electron charge, respectively. Others 
are material-dependent parameters, e is the dielectric constant. Vb is the built-in potential. 
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Na and A^d are the density of acceptors and donors, respectively. We use Na = Nb = 
where n 2 D is the 2D electron density on the Te-top domain. n 2 D is estimated by the Fermi 
wavenumber and the wavenumber of the electron standing wave at 0 mV, gin. Namely, 
^20 = k^^/{2 t:) and fcp = gm/2. Since gin ~ 2.8 nm“^, we get Na = iVo 1-7 x 10^° cm“^. 
Other parameters are e ~ 19 (Ref. 42) and Vb 1.6 V that is estimated from the difference 
of the work function shown in Fig. 4(e). By using these values, we get d ~ 6.2 nm. 

Appendix F: Apparent height in a topographic image 

Here we give a simple mathematical description about multiple contributions to apparent 
height of a constant-current topographic image. The STM tunneling current / is given by 

feV 

I(r,z,V) = CeM-2K(r)z) N(r,E)dE (FI) 

Jo 

while z is the tip’s surface-normal coordinate, V is the sample bias, and N{r, E) is the 
sample’s LDOS at lateral locations r and energy E. k is related to the tunneling work 
function nir) = ^J2m(f){r)/h. e, m, and h are the elementary charge, the electron mass, and 
the Planck constant, respectively. Given this formula, a constant-current (Jq) topographic 
image Zcdr; Iq, Vq) af a sample bias voltage Vq is given by 

h (C \ 

zccir; /o, W) = , - In - / N{r, E)dE . (F2) 

,j8m(p{r) Vo Jo J 

Eq. (F2) indicates that apparent height of a constant-current topographic image is low when 
the work function is large or the LDOS is small. This is a consequence of a feedback loop 
that gets a tip closer to a sample to compensate fast decay of the wave function or reduction 
of LDOS available for tunneling. In the case of BiTel, the work function is larger in the 
I-top domain as shown in Fig. 4(e). The LDOS is likely to be smaller in the I-top domain 
because the 2D hole state in the I-top domain is caused by strong inversion whereas the 2D 
electron state in the Te-top domain by accumulation. Namely, both the work function and 
the LDOS contribute to the apparent height lower in the I-top domain. 

Appendix G: No extrinsic surface modification 

Here we describe that our samples have neither surface absorbates nor additional defects 
possibly changing carrier density at the surface. As is evident in the topographic images 
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shown in Fig. 1, no surface absorbate exists. Surface defects additionally generated by 
cleavage are excluded by counting the density of defects as follows. 

To evaluate the density of defects near the surface, we assume that observed defects he in 
the top unit layer. This assumption gives the highest estimate of defect density. By counting 
defects in several topographic images, we estimate that the density of these defects is about 
5 —8 X 10^^ cm“^. This value, slightly larger than the bulk electron density (2.5 x 10^^ cm“^), 
is just as expected because a considerable amount of acceptors exists in this material as well 
as predominant donors^^. Therefore, the density of defects near the surface is virtually 
identical to that in the bulk without increased by cleavage. 
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